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ABSTRACT 


This study investigates biodiversity and structure of spring diatoms associations from Sebou Basin, one of 
Morocco's largest superficial hydrosystems. The physico-chemical results reflect both anthropogenic and 
natural pollution, mainly from mineral sources. 226 taxa were identified, 92% are pennates. The result of 
Shannon and equitability indices showed that diatom populations were generally diversified, which depends on 
degree of specific richness at each station. Thus, the higher number of taxa in site Og (120 species) reflects an 
ecologically more stable and more diversified ecosystem dominated by E. minuta, opposed to S, (Sebou) and 
B» (Beht), which have an unbalanced community represented by smaller numbers of taxa (48 in S4 and 52 in B;) 
strongly dominated by N. capitellata and N. palea. The Canonical correspondence analysis distinguished three 
groups: group 1 with SO, B1, Og and In located far from pollution sources, group 2 composed of polluted sites 
(S1 and B2), while group 3 with of S2, S3, and S6 are located in the lower section of the Sebou Basin. Multi- 
response permutation procedures confirmed the three site groupings. It was found that the water quality located 
downstream from the urban rejections is degraded, indicator species analyses showed a succession of changes 
in diatom species from upstream sites to downstream. Six species were determined as indicator species; 
Groups 1 and 2 only had one species each, G. olivaceoides (8696) and H. amphioxys (8796), respectively, while 
group 3 had four species, N. recta (93%), N. sigma (95%), N. tripunctata (96%) and D. oblongella (100%). 


Keywords: watershed of Sebou (Morocco); diatoms; multi-response permutation procedures (MRPP); indicator 


species analysis (IS); Canonical Correspondence Analysis (CCA). 


INTRODUCTION 


Diatoms colonize a wide variety of habitats and are an 
important part of the biomass in most aquatic 
ecosystems. They have long been used as the best 
biological indicators of Europeans waters quality (Prygiel 
and Coste 1999; Yallop and Kelly 2006; Rimet 2012). 
Thus, precise knowledge of ecological requirements of 
certain taxa helps the development of simple indicators 
systems that characterize and quantify water pollution. in 
contrast, morocco's diatoms were many studied, but 
among aquatic phytoplankton and/or in limited areas 
(Fqih Berrada, Berrada, and Benzekri 1999; El Ouahabi 
et al 2007; Siddour et al 2007, Loumrhari et al 2009; 
Gallouli et al 2014). Moreover, morocco's diatoms 
ecological studies refer to European ecology, while their 
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exact range of ecological conditions of where they occur 
in Morocco are still unclear. The most important works 
performed in Morocco are those made by Fekhaoui, 
Hamada, and Dakki (1988) on Sebou River; on Oued 
Tensift and its tributaries (Maiffi 1988; Cazaubon and 
Badri 1994): and on wadis Hassar and Mellah rivers 
(Fawzi 2002). 
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Figure 1: A. Geographical situation of the Sebou River watershed of (MARA 1970). B. Map showing the location 


of the study sites 


However, no studies have been published for the 
diatoms found in Sebou Basin that encompasses one of 
largest superficial waters networks and that suffers 
sustainable socio-economic development, leading to 
continue exposure to various pollution (agricultural, 
industrial and urban) and prejudicing on aquatic life. The 
objects of this study are: 

- Analyze of diatoms and physico-chemical content 
from this great Moroccan river and its most important 
tributaries (Innaouéne, Ouergha and Beht). 

- Then estimate of ecological preferences of abundant 
diatoms in order to enrich the bases of Moroccan given. 


MATERIALS AND METHODS 


Study sites: 

The Sebou Basin (Figure 1A) occup an area of about 
40,000 km”, it is located between 33° and 35° North, and 
4? and 7? West. Sebou River, is a principal water course, 
originates in the Middle Atlas and flows mainly across 
Lias carbonated formations (Fedan 1988; Charriérer 
1990; Laadila 1996; Dridri and Fedan 2001a, 2001b). 
Rare Paleozoic buttonholes and Quaternary basalts 
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outcrops interrupt the area's monotony. The Gharb plain 
contains thick fluvio-marine beds and sandy-clay 
alluviums deposited during the Quaternary (Rif, Middle 
Atlas and Prerif) (Azzaoui 1999). 

Otherwise, roughly a fifth of Morocco's total 
population, 2296 of the national rural population and 1796 
of the urban population live in the Sebou Basin. Similarly, 
it is important region of the country for its large 
agricultural surface area, and a diverse industrial 
complex: oil, sugar and paper mills, tanneries, textiles, oil 
refining, etc. 

In addition, the Sebou Basin contains nearly a third of 
the nation's surface water (33%) and 20% for ground 
water (ABHS-M 2006). 

The studied stations were selected taking into account 
the: 

- Relationship to pollution sources: samples collected 
upstream and downstream of discharges and then some 
stations were considered reference stations because of 
the location far from pollution discharges. 

- Ease of access and sampling. A set of 11 stations, 
located according the longitudinal gradient of the river 
course (from SO to S6) (see Figure-1). 
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- S0: Sebou River, about 20 km upstream from the 
Fez city. No notable anthropologic activity (troughs or 
washing). 

- S1: Sebou River, to 1.5 km downstream of the Oued 
Sebou confluence with Oued Fes. Downstream from 
the city of Fez. 

- |n: Innaouéne river (tributary of Wadi Sebou), about 5 
km downstream of ldriss 1st water-dam. This is a 
reference station. 

- $2: Sebou River, about 8 km upstream from the town 
of Kariat Ba Mohammed. 

- Og: Ouergha River (great Oued Sebou tributary), to 
10 km upstream of Khnichat and about 20 km 
downstream of the dam El Wahda. This is a reference 
station. 

- S3: Sebou River, to 3 km downstream from the city 
Mchraa Bel Ksiri. 

- $4: Sebou river; just downstream of the town of Sidi 
Allal Tazi. 

- $5: Sebou River, just downstream retained of the 
dam custody and influenced by the tides. 

- B1: Beht River (largest tributary of Wadi Sebu), to 10 
km upstream of the releases of the town of Sidi 
Slimane. No of the anthropologic activities except 
washing and watering by villagers. 

- B2: Beht River, located just 
discharges of domestic wastewater. 

- $6: Sebou River estuary, upstream of the mouth and 
downstream of the domistic and industrial wastewater 
of city of Kenitra. 


downstream of 


Physical and chemical analysis: 

Five parameters were measured jin situ: water 
temperature, pH and conductivity with a WTW apparatus 
while salinity and dissolved oxygen were measured with 
a salininometer and oximeter, reference 320. The 21 
water samples were collected in polyethylene bottles and 
stored in an icebox (4?C) until arrival at the laboratory. 

Fourteen physic-chemical parameters were measured 
and dosed. Chemical oxygen demand (COD) was 
analyzed immediately (within 24 hours of sampling), 
while biochemical oxygen. demand (BOD;) was analyzed 
followng AFNOR standards. The other parameters 
(alkalinity, nitrite, nitrate, ammonium, ortho-phosphates, 
sulfates, chlorides, total hardness, calcium, magnesium, 
sodium and potassium) were analyzed respecting their 
shelf life using methods prescribed by AFNOR standards 
or Rodier et al. (2009). 


Diatoms sampling, treatment and determination: 

Diatom samples were mainly epilithic and were 
collected from the different sampling stations (Table-1), 
according to the standard NF T90-354 (AFNOR 2000, 
2007). The top of several submerged stones were 
scraped (approx. 200 cm?). The samples were 
immediately fixed with 10 ml of 1096 formalin. In the 
laboratory, diatom samples were boiled in H20; and then 
in HCl, before thoroughly rinsed. Permanent slides are 
produced by mounting a few drops of the prepared 
suspension with a drop of Canadian balsam. At least 400 
individuals were counted at a magnification of 1000 x 
using oil immersion. The systematic identification of 
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diatoms was performed based on the taxa key found in 
the standard NF T90-354 AFNOR (2003, 2007); 
Bourrelly (1968); Compére (2000); Lange-Bertalot (2000) 
and DRIEE-IF (2013). 


Diversity and equitability Indices: 

The Shannon-Wiener diversity index (H') expressed in 
"bit", is amount of information provided by a sample on 
the population and on the species structures (Daget 
1976).It is calculated using the following formula: 

'z -5 (ni/ N) * Log 2 (ni/ N) 
H': species diversity. 
N: total populations of species. 
ni: Population size of species i. 
Equitability also called relative or regularity diversity 
compares the stand structures. It is defined as the ratio 
of real diversity to the maximum diversity (H max) and is 
expressed as follows: 
E- H'/ Hmax with Hmax- Log2 (S) 
Where S is the number of species forming the stand. 

Remember that Shannon-Wiener index (H') varies 
between 0 and H max (Barbault 1995). A diverse 
population should have elevated species richness and 
equal distribution. 


Statistical analysis: 

Canonical Correspondence Analysis (CCA), multi- 
response permutation procedures (MRPP) and indicator 
species analysis (IS) were carried out to assess diatom 
species-environmental relationships (Borcard, Legendre, 
and Drapeau 1992; Bere and Tundisi 2011b; Birks 
2012b; Bertrand, Serieyssol, and Ector 2015) using taxa 
having 2196 of the community and occurring in more than 
one sample. The following physico-chemicals were used 
in the analysis: dissolved oxygen (DO), conductivity 
(Cond), sulfate (SO4), phosphate (PO,), nitrate (NOs), 
ammonium (NH4) and pH. However, conductivity data 
was log 10 transformed prior to analyses. The statistical 
methods were performed using PCORD 6 (McCune and 
Mefford 2006) and the following CCA options were 
selected: the options were columns score standardized 
by centering and normalizing, scaling of ordination — 
comprise representation of species and samples, scores 
for graphing samples scores are linear combination of 
factors, randomization test on, test null hypothesis — no 
relationship between matrices and number of runs 999. 
Data from stations S4 and Ss were not exploited for their 
physical and hydraulic characteristics because of poor 
diatom colonization. 


RESULTS AND DISCUSSION 


Physical and chemical variables (Figure 2, 3, 4, 5 
and Table 1): 

The results of all physico-chemical variables 
measured at 11 sites are given in Table 2 and Fig 2; 3; 4; 
and 5. The changes within basin are generally related to 
nature of substrate but are also related to the urban, 
industrial and/or agricultural situation. Water temperature 
varied between 15.5°C at site B1 to 22.7?C at S6 site. 
The pH varied entre 7.38 at B2 to 8.9 at S4, indicating 
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Figure-2. Canonical correspondence analysis (CCA) performed on the diatoms collected in nine stations of 
the Sebou Basin related to major variables. Cond: conductivity; O2: dissolved oxygen; PO: 
orthophosphates; NO3-: nitrates; NH4+: ammonium; HCO3-: bicarbonate and pH. 


Axis 2 


Axis 1 


alkaline waters. It is due mainly to calcareous marl and 
limestone bedrock (Fekhaoui and Pattee 1993; EI Blidi 
and Fekhaoui 2003), combined of  insolation 
(temperature), and increase of plant photosynthetic 
activity in spring (Fawzi 2002; Benabdellouahad 2006). 
The conductivity varied between 533 us/cm at Og to 
1402 us/cm at S1, respectively, indicating a moderately 
accentuated mineralization to elevated (Rodier et al. 
2009) The maximum values of the majority of 
parameters such as COD (113mg/L); BOD (130 mg/L); 
HCO; (415 mg/L); CaCO; (360 mg/L); Ca” (132.26 
mg/L); Mg* (227.74 mg/L) K* (15.60 mg/L); PO,” 
(164.42 ug/L); and NH,” (192.12 mg/L) were observed in 
site S1 with minimum oxygen-content of 0.29 mg/LO;. 
The sites S0; In; Og; and B1 were the most oxygenated, 
and less rich in organic and, or mineral parameters. High 
nitrate values were measured in station S2 (5.67 mg/L of 
N), with a nitrite peak observed at station S3 (0.86 mg/L 
of N). However, conductivity and mineral ions analyses 
(Cl; SO4^; Na*; CaCOs; HCO;; Mg*'; Ca?) also indicate 
a lithological origin. 

The adjacent agricultural lands and the proximity to 
urban centers of some stations amplify mainly the 
sodium and sulfate ions, and also leads to increased 
chlorides. Makhoukh et al. (2011) and Libiad, Khabbach, 
and Ennabili (2012) made the same remarks in their 
respective studies of the Moulouya and Innaouéne 
hydro-systems. 

The statistical study showed that the first three axes 
of the CCA (Fig 2) explained 59.3 96 of the data (axis 1 
22.896, axis 2 18.5%, axis 3 18.0%). Three main groups 
were found; group 1 with SO, B1, Og and In, located far 
from pollution sources (reference sites). In contrast, 
Group 2 composed of polluted sites, S1 and B2, are 
located downstream from Fes and Sidi Slimane cities 
respectively, while group 3 with of S2, S3, and S6 are 
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located in lower section of Sebou Basin. For group 2, 
samples are situated in the left half of the graph (Fig 2) 
indicated by high conductivity and high amounts of 
carbonate (HCO3), ammonium (NH4) and phosphates 
(PO). The upper right quadrant had group 1 samples with 
high dissolved oxygen values while group 3 samples are 
situated in the lower left quadrant with high amounts of 
nitrates. 

The graph of the first two axis scores (Fig 3) show 
some of the relative changes that occur from headwater 
to the ocean. As you move down river starting with group 
1 having low conductivity, low concentrations of HCO;, 
PO, NH, and NO; with high levels of dissolved oxygen to 
Group 2 has high conductivity and high levels 
concentrations of HCO3, PO and NH,, increasing levels 
of NO; and lower levels of Op and pH, followed by group 
3 with a return to lower conductivity and concentrations 
of HCO3, PO and NH4, but with high concentration of 
NO3. 

The chemical variation within the three groups (1- 
upstream, 2-polluted sites and 3-downstream) are 
shown on Figure 4. Groups 1 and 3 have similar profiles 
with S6 showing increased salinity caused by tidal 
fluctuation. Within group 2, site S1 showed the greatest 
amount of pollution. 

In Figure 5, the sites are position in the geographic 
distance from upstream to the sea and variation of the 
organic pollution and chemical variables are observed 
along with the different CCA groups the sites belong to; 
group 1 — upstream, group 2 — most polluted and group 3 
— downstream. 

Interestingly, when the upstream tributaries (In and 
Og - group 1) enter the Sebou River, there is a lowering 
of all of the variables. However, BOD, PO, sulphates 
(S04) and chlorides (Cl) increase at S1 after Oued 
Innaouéne (In) waters enter the river but when Oued 
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Figure-3. Canonical correspondence analysis (CCA) axis scores showing the changes in ecological 
condition occurring between groups 1, 2 and 3 or from upstream to polluted area and then to downstream 


zone of the Sebou Basin. 
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Figure-4: Variation of organic pollution and nitrogenous variables belonging to the different groups (1 — 
upstream, 2 — pollution and 3 — downstream). Dissolved oxygen (OD mg.L-1), chemical oxygen demand 
(COD mg.L-1), biological oxygen demand (BOD mg.L-1) HCO3 (mg.L-1) Ca CO3 (mg L-1) Ca (mg.L-1), Mg 
(mg.L-1), Cond (ps.cm-1), CI (mg.L-1), Na (mg L-1), K (mg L-1), SO2 (mg L-1), PO (pg.L-1), NO3 (mg.L-1), 


NO2 (mg.L-1), and NH4 (mg.L-1). 
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Ourgha water (Og) enter the river, all of the variable 
seem to decrease from the upstream site (S2) to the 
downstream site (S3). Interestingly, even the relatively 
high BOD value had not effect on site S3. There is 
however an increase in Conductivity, Cl, Na, CaCOs, 
Calcium (Ca), HCO3 at the estuary site (86). Chemical 
trends for the Oued Beht tributary could not be 
determined considering the large distance between the 
two upstream sites (B1 and B2) and the estuary site 
(S6). 

Besides that, the maximum values of mineral and 
organic elements were recorded in Group 2 sites (S1 and 
B2) reflecting the large amounts of pollution discharged 
upstream from these sampling stations. This caused 
excessive oxygen consumption creating critical de- 
oxygenation, accompanied by a pH drop. The reduction 
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of organic matter accompanied by re-oxygenation of 
waters in group 3, show river self-purification and the 
dilution effect of well-oxygenate waters entering the river 
coming from cleaner tributaries (In; Og and B1- group 1). 
Behavior of nitrogen compounds (NO3, NOs, and NH) in 
the basin is related to the reaction of oxidation- 
reductions. Indeed, the nitrate increase is attributed to 
nitrification reactions of ammonia nitrogen via nitrite 
nitrogen (Jaghror 2013). Conversely, a decrease 
indicates ammonification. Note that the domestic and 
industrial wastes from Kenitra city are less noticeable in 
the Sebou Estuary (S6-group 3) in comparison with Fés 
(S1) or Sidi-Slimane (B2). Indeed, the waters of the 
station S6 are subjected a continuous stirring by tides, 
which disperses the pollution and/or is absorbed by soil. 
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Figure-5. Spatial variation of organic pollution and nitrogenous variables along with the different groups 
(1 - upstream, 2 - pollution, and 3 - downstream) A. The Sedou River and two tributaries (O. Innaouéne 
and O. Ourgha), and B. the O. Beht tributary which enter the Sebou River near site S6. The sites are placed 
in the order as they occur along the river from upstream to downstream. Dissolved oxygen (OD mg.L-1), 
chemical oxygen demand (COD mg.L-1), biological oxygen demand (BOD mg.L-1), SO2 (mg.L-1), PO 
(Mg.L-1), NO3 (mg.L-1), NO2 (mg.L-1), NH4 (mg.L-1), CI (mg.L-1), Na (mg.L-1), K (mg.L-1), HCO3 (mg.L-1), 
Ca CO3 (mg.L-1), Ca (mg.L-1), Mg (mg.L-1), and Cond (ps.cm-1). 
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Diatoms analysis 


Specific richness, index of diversity and 
evenness (Table 2 and 3): 

The spring diatom flora of Sebou Watershed is rich 
and diverse. 226 species belonging to 57 genera and 5 
orders were identified in nine samples. The 
Pennatophycidae had most represented (209 species 
classified in 3 orders and 50 genera) while the 
Centrophycidae had only 17 taxa belonging to 7 genera 
(only 7% of the harvested diatom). Nitzschia followed by 
Navicula (47 and 31 species, respectively) had highest 
species diversity, with 36% of identified species. 
Gomphonema and Surirella had respectively 17 and 11 
species; the other 53 genera had fewer species with 22 
having only one. The most common species throughout 
the Sebou Basin were: Nitzschia palea (Kütz.)W.Sm., 
Nitzschia frustulum (Kütz.) Grunow var. frustulum and 
Nitzschia dissipata (Kütz.) Grunow var. dissipata were 
observed in all sites. The population with highest diatom 
richness was found at site Og with 120 species (Table3 
and 4), that presents 52.496 of the diatom community 
collected in the Sebou Basin. To the contrary, stations 
S1 (Oued Sebou) and B2 (Oued Beht) had the lowest 
number of species (48 and 52 respectively). 

Twenty-eight species had a relative abundance 2 5% 
(table 3) while 92 species had a relative abundance of 
1% in at least one sample, indicating that the 
environment was diverse. The results of Shannon- 
Wiener index (H') and equitability (E) have same species 
richness trends for all the stations without showing a high 
analogy between stations S3 and B1, and do not exhibit 
significant spatial differences (with the exception of S1 
and B2), indicating that the environment is generally 
diversified. The diversity degree depends on the species 
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richness of each station. The obtained values of diversity 
ranged between 2.861 for S1 and 5.732 for Og (Table 4), 
the evenness varied from 0.512 at S1 to 0.840 at S3 and 
B1. Thus, the highest richness in Og station reflects an 
ecologically stable and diverse aquatic ecosystem 
including 73 taxa having lower relative abundance than 
196. At SO and B1, H' is close to H' max with an diatom 
community less diverse and less stable than Og. This 
difference is most likely related to the water depth and 
high current velocity, as all these stations are located 
away from human disturbances. Indeed, according to our 
field observations, the water depth at Og was deeper in 
comparison with SO, In and B1 sites. The sites S1 and 
B2 (Group 2) showed the lowest diversity indices values 
revealing the presence of unbalanced communities and 
are represented by a small number of taxa. These 
diversity index values confirm deterioration of water 
quality at the two stations, which are dominated by 
Nitzschia capitellata Hust. described in Europe as 
indicator species of polysaprobic waters, hypereutrophic 
(Krammer and lange- Bertalot 1988; Van Dam, Mertens, 
and Sinkeldam 1994; Hofmann, Werum, and Lange- 
Bertalot 2011), and rich in electrolytes (Taylor, Harding, 
and Archibald 2007), which is the case for these two 
sites. Fawzi (2002) also noted in his study a major 
development of Nitzschia capitellata in moderately 
alkaline waters, eutrophic, and heavily polluted by 
ammonium. The increase of both diversity indices and 
species richness in S2 and S3 (Group 3) indicates that 
the population resumes a balance following water quality 
improvement by the self-purifying capacity of the river. 
The decrease of diversity index and number of taxa 
colonizing near the estuary (site S6) are due to the 
reappearance of taxa living in eutrophic conditions 
prefering brackish waters to salt waters such as 
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Nitzschia filiformis (W.M. Smith) Van Heurck var. 
filiformis (Van Dam, Mertens, and Sinkeldam 1994), 
Cyclotella  meneghiniana  Kütz (Hürlimann and 
Niederhauser 2007; Taylor, Hading, and Archibald 2007; 
Houk, Klee, and Tanaka 2010; Bey and Ector 2013) and 
Nitzschia clausii Hantzsch (Germain 1981; Van Dam, 
Mertens, and Sinkeldam 1994; Bey and Ector 2013). The 
S6 site periodically receives the rising tide from 
downstream and the upstream waste rejects. However, 
the real degree of chemical pollution, essentially salinity, 
compared with organic pollution cannot be full 
determined as samplings occurred at low tide. The 
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community of group 3 (S2, S3 and S6) indicated 
freshwater to slightly brackish, even marines conditions, 
oligosaprobic to  mesosaprobic, eutrophic, rich in 
calcareous and electrolytes as underlined by the authors 
named above. The main taxa were Craticula halophila 
(Grunow ex Van Heurck) Mann (Krammer and Lange- 
Bertalot 1986; Taylor, Harding, and Archibald 2007), 
Navicula reichardtiana Lange-Bert. var. reichardtiana 
(Krammer and Lange-Bertalot 1991b; Lange-Bertalot 
2001; Taylor, Hading, and Archibald 2007; Lange- 
Bertalot et al. 2011), Navicula menisculus Schumann 
var. menisculus (Lange-Bertalot 2001), Nitzschia sigma 
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(Kütz.) W.M.Smith (Ziemann 1991; Busse, Jahn, and 
Schulz 1999; Fawzi 2002) and Navicula gregaria Donkin 
(Kelly 1998; Lange-Bertalot 2001; Fawzi 2002). 


Statistical analysis (Fig 6): 

The Canonical Correspondence Analysis showed 
that taxa preferred water with high dissolved oxygen 
concentrations were Achnanthidium  jackii, Amphora 
pediculus, Cyclotella ocellata Pantocsek, 
Gomphosphenia lingulatiformis, Navicula menisculus and 
Surirella helvetica. Cocconeis pediculus was associated 
with chloride concentrations and Cyclotella meneghiana, 
Navicula gregaria and Nitzschia filiformis was observed 
in areas with nitrates concentrations. Water with high 
concentrations in phosphorous, ammonium, bicarbonate 
and conductivity were linked with the following species: 
Nitzschia capitellata, Nitzschia palea, and Tryblionella 
constricta (Kütz.) Poulin. Nitzschia Palea is known to be 
associated with higher concentrations of phosphorus and 
nitrogen (Chaïb et al, 2011). 

Using the three groups determined by the CCA 
analysis,  Multi-Response Permutation Procedures 
(MRPP) was performed to verify the position of Og in 
group 1. The greater the T value the greater difference 
between the groups. The p values indicate that groups 
are statistically valid (Table 4). MRPP was run three time 
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with site Og being added to one of three groups. A T 
value of -4.32 was determined when Og was associated 
with group 1 (upstream group) while when placed within 
group 3 (downstream group) a lower T value of -4.09 
was observed. When associated with group 2 (polluted 
group) non-significant result was obtained (Table-4). 
Indicator Species analysis (IS) was run using the 
three groups. It should be noted that those species found 
in almost all of the samples were excluded from the 
analysis. Result found that certain species occurred in 
certain groups in greater than 80%. This marks species 
composition changes that occurred as moving 
downstream from group 1 to group 3 (Table-5). 
Interestingly, five species occurred in all three analyses, 
relative abundance, relative frequency and indicator 
species. In the relative abundance (percent of each 
species in each group defined by topographic position), 
seven species (Achnanthes coarcta (Breb) Grunow, 
Fragilaria capucina Desmazieres var. capucina, 
Gomphonoma clavatum Ehrenb., Gomphonema 
angustatum (Kütz) Rabenh., Surirella biseriata Breb. in 
Breb. and Godey, Gomphonoma  olivaceoides Hustedt, 
and Gomphonoma gracile Ehrenb.) were found in > 80% 
in group 1, three species (Gomphonema angustatum, 
Surirella biseriata, and Gomphonema olivaceoides) were 
present in group 2 and one, Gomphonema gracile, in 


Figure 6: Canonical correspondence analysis (CCA) axis scores showing the changes of community 
occurring between groups 1 (dotted circle), 2 (black circle) and 3 (grey circle) or from upstream to polluted 
area and then to downstream zone of the Sebou Basin (code names are in Table 3). 
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Table-1: The physicochemical parameters analyzed in the study sites 


Station] 0D | COD | BOD | HCOS | caco: | Ca^ | Mg” | Cond CF | Na | K ] SOx] POS | NO: | NO; | NH" | cc | py 
(mg/) | (mg/l) | (mg/l) | (mg/l) (mg/l) (mg/l) (mg/l) us/Cm | (mg/l) | (mg/l) | (mg/l) | (mg/l) | (ug/l) (mg/l) | (mg/l) | (mg/l) 

so | 10.15 | 40 0 255 280 | 11222 | 16778 | 994 355 | 138 | 1170 | 25.95 | 1988 | 198 | 003 | 000 | 203 | 81 
sı | 029 | 113 | 130 | 415 360 | 13226 | 22774 | 1402 | 5325 | 138 | 1560 | 2579 | 16442 | 151 | 021 | 92! | 20 | 76 
In. | 1067 |. 17 0 171,86 149,1 zo| 295 | 4723 | 079 T 21 | 848 
$2 | 759 | © 30 2 | 17980 | 752 | 5609 S0 | 2654 | 8875 | 567 f ; 178 | 83 
Og | 1036 | 70 0 | 107,82 | 533 142 60 | 1404/ 1459 | 062 191 | 817 
s+ | 98 | 0 25 19279 | 581 | 1704 30 |1154| 6415 | 184 168 | 88 
sa | 914 | 7 15 123,85 560,9 70 | 81} 2594 | 196 166 | 89 
s | 938 | 24 | 25 | 20784 | 588 355 80 |2651| 3048 | 084 165 | & 
B1 | 14.97 | 45 0 16004 | 815 | 7455 904 | 16,10 | 256 i 155 | 7,61 
B2 | 824 | 54 | 50 17427 | 880 | 1917 60 | 43.17 | 5848 | 252 167 | 7.38 
se | 447 | 21 | 2222 175.79 3763 | 161 | 390 | 8373 5.22 227 


Table-2. List of abondantes spring diatoms collected with relative abundance (%) in sites study of the Sebou 
Watershed. 


| Abbrevation | So | Sı | In | & | Og | ss | B | Be | Se | 
Pennates 


Ee E 
A Jacki Rabenh. ær |18|03|23]| 0 |33|23 | $5 [05 | 9 
A. minutissimum Kitz. Czarnecki | ADM | 85] 05 | 0 | 03 | 23 | 13 | 4 | 02 | 05 | 
Anphore iya Eh. aus —05|03]| 0 | 03. | 03] 0 | 5 [05 | o 
A pedes Kt) Gruno [ aep | 3 [0s[0s|03 | 18| 1 | t5 [05 | o 
| C. pediculus Ehrenb_ | PED [13] 1 | 0]|05]|25]|15| 4 | 151] 
ER umi A ne SR NU M CN 
| C. halophila (Grunow ex Van Heurck) Mann  . | | | CHAL |03| 0 | 0 | 03 | 03 | 75 | 0 | 0 | 
UL TOI net Our pipe is oa [ét FL ERES 
nn (esr Kranmervar Te | "CH | o UE ee a 


Encore minuta (Hilse in Rabh.) D.G. Mann in Round Crawford and ENMI eA el eel ee | 


| Encyonopsis microcephala (Grunow) Krammer 1997 — | ENCM [03,03 | 258] 08 | 35 | 15 | 0 | 02] 0 | 
| G. intermdiumHustedt | aM [os| o| of 1]|03]| 0 | 78 | 12] 0 | 
| G. pumilum (Grunow) ReichardtandLange-Bet — 1 | GPUM |03| 0 [138|05]| 1 | 0 | 33] 1 | 03 | 
E Gomphosphenia Ingulatformis (Lange Bent. and Reichard) Lange Ben | GPi | os| o | 0 | i |025| o |675 | 125 | 0. 
| N. gregaria Donkin | NGRE | 005] 0 | 1 | 13 |123| 05 | 42 [118] 
| N. menisculus Schumann var. menisculus "| NMEN [03[|03][05/|68]| 1] 1 [03] 0 | 03| 
| N. reichardtiana Lange-Bert. var. reichardtiana "|  NRCH | 0 | 0 | 03 | 18 | 0.8 | 13 | 0 | 0 | 0 | 
| N. amphibia Grunow f. amphibia "D! NAMP _]01]01|]11}23]03] 0 | 23 [99 | 03 | 
LN. capitellata Hust. in A Schmidt and al." |  NCPL | 0 |563] 03 | 0 | 03] 0 | 0 [387] 03 | 
| N. clausii Hantzsch | NCLA | 0 | 0 | 0 | 0 | 13] 0 | 0 | 0 | 58] 
LA filformis (W.M Sri) Van Heurck var fifommis | NL | 0 [93| 9 | 0 [03| 4 | 0 | 9 [ 14, 
| N. frustulum(Kutzing) Grunowvar. frustulum "| NFR [28] 25] 1 | 85 | 45] 58 | 53 | 37 | 13 | 
|Npalea(Kitz)WSm "| NPA | 15] 98 | 08 | 15 | 08 | 48 | 35 | 129] 08 | 
| N. thermaloides Hustedt "| NTE | 0 | 85] 0 | 0 | 0 | 0 | 0 | 0 | 0 | 
| Surirella angusta Kütz bo | SANG [13| 0 [| 0 [08] 03] 0 | 6 | 0 | 55] 
[S antoquiensis Sala, Ramirez, Plata-Diaz andVVoulloud | suan [ss] o | 0 | 0 | 0 | 08| 0 [92[ o] 
| S. helvetica Brun  — O O 1 | 1 O SHE [333] 05 | 03 | 35 | 28 | 05 | 75 | 0 | 0 | 
| T. constricta (Kütz.) Poulin ^ ^ ^ à à à | mco | 0 | 03 | 0 | 08 | 13 | 13 | 03 | 79 | 0 | 
CE ee ie 
| Conticribra weissflogii (Grunow) Stachura-Suchoples and Williams — | =CTWE | 0 | 0 | 0 | 0 | 0 | 0 | 63 | 0 | 0 | 


[CmeneghineraKitz L| CMN | 0 fos fos os [93 | t | o | * | 5] 
| C. ocellata Pantocsek — 1 1 1 | | [| COCE | 0 | 03 | 10 | 38 | 48 | 3 [33] 0: |. 25 | 


group 3. As for relative frequency (the percent of 
samples in each group in which a species is present), 
five species (Gomphonema  olivaceoides, Navicula 
cryptotenella Lange-Bert., Gomphomema parvulum Kütz 
var. parvulum f. parvulum, Amphora pediculus (Kütz.) 
Grunow, and Cymbella affinis Kütz .var. affinis) were 
found to occur at 100% in group 1 but two of the species 
(Gomphonema parvulum var. parvulum f. parvulum and 
Amphora pediculus) also occurred at 100% in group 2 
along with Hantzschia amphioxys, Nitzschia fonticola 
Grunow in Van Heurck, Nitzschia dubia W.M.Smith and 
Craticula accomoda (Hustedt) Mann at 10096. Group 3 
had nine species (Cymbella affinis Kütz.var. affinis, 
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Table-3: Diatom diversity indices for the sampling sites: Shannon index (H’), index of maximum Shannon 


(H'max), and the equitability (E) 


S0 $1 In $2 Og $3 B1 B2 S6 

Total abundance 62 48 70 120 82 58 52 67 
H' 4.220 2.861 4.195 4.999 5.732 5.338 4.922 3.579 4.914 
H'max 5.977 5.585 6.129 6.340 6.907 6.358 5.858 5.700 6.066 
E 0.706 0.512 0.684 0.789 0.830 0.840 0.840 0.628 0.810 


Table-4. Multi-response Permutation Procedures results. A. Og belonging to group 1 (upstream sites) with a 
total analysis having a T value z -4.32. B. Og belonging to group 3 (downstream sites) with total analysis 
having T value = -4.09. C. Og belonging to group2 (polluted sites) had non-significant results. 


A Group 1 | Group 2 | Group 3 B Group 1 | Group 2 | Group 3 
Group 1 -2.86 -2.58 o Group 1 -2.2 2.32 o 
Group2 | 0.0176 -2.23 a Group 2 | 0.0001 -2.99 S 
Group 3 | 0.0145 | 0.0001 ES Group3 | 0.0221 0.0157 ü 
p value 


Mayamaea atomus (Kütz.) Lange-Bert., Nitzschia recta, 
Nitzschia sigma, Navicula tripunctata, Actinocyclus 
normanii (Greg. ex Grev.) Hust. morphotype normanii, 
Navicula veneta Kütz., Cyclotella meneghiniana Kütz. 
and  Diploneis oblongella (Naegeli) Cleve-Euler). 
However all but Diploneis oblongella (Naegeli) Cleve- 
Euler were also found in group 1 and Cymbella affinis 
Kütz .var. affinis also occurred in group 1 at 100%. Three 
species not only occurred with relative frequency group 
1, but also in group 2 (A. normanii, N. veneta and C. 
meneghiniana). However, Diploneis oblongella occur 
only in group 3 at 100%. Finally, only six species were 
found as indicator species (percent of perfect indication 
of each species for each group). 

Groups 1 only had one species, Gomphonema 
olivaceum var. olivaceoides (Hustedt) Lange-Bertalot 
with 86% (Table 5) which preferred oxygenated sites, 
poor in organic and mineral matter, the same ecology 
was defined by Bey and Ector (2013) in their study of the 
rivers from the Rhone-Alpes region. The analysis of 
diatom assemblages underlined that taxa present in 
group 1 are generally cosmopolitans, colonizing 
freshwaters to slightly brackish, oligosaprobic, oligo- to 
mesotrophic but preferring the calcareous waters with a 
moderate to high electrolyte contents, according to Van 
Dam, Mertens, and Sinkeldam 1994; Taylor, Harding, 
and Archibald (2007); Hofmann, Werum, and Lange- 
Bertalot (2011); Bey and Ector (2013). The only indicator 
species of group 2 is Hantzschia amphioxys (Ehr.) 
Grunow described by Van Dam, Mertens, and Sinkeldam 
(1994); Hofmann, Werum, and Lange-Bertalot (2011); 
and Bey and Ector (2013), as an indicator of polluted 
environments. Furthermore, the main diatoms of group 2 
were representatives of the genera of Nitzschia, 
indicators of heavy metals contamination (Guasch et al. 
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2009; Tlili et al. 2011), salinity (Hofmann 1997; Rimet 
2009), eutrophic (Rimet, Berthon, and Bouchez 2010) 
and / or organic pollution (Kobayasi and Mayama 1982; 
Van Dam, Mertens, and Sinkeldam 1994; Martin et al. 
2010; Trobajo et al. 2013). while group 3 had four 
species, Nitzschia recta (93%), Nitzschia sigma (95%), 
Navicula tripunctata (96%) and Diploneis oblongella 
(100%). therefore are characteristic of polluted sites by 
nitrates with high mineralization. Note that among the six 
indicator species determined, as we consider Diploneis 
oblongella as the best indicator of group 3 since it was 
present with a frequency of 100%. Four taxa in this group 
had a indicator value superior to 80% were Nitzschia 
recta Hantzsch in Rabenhorst; Nitzschia sigma (Kutzing) 
W.M.Smith; Navicula tripunctata (O.F.M_ler) Bory. and 
Diploneis oblongella (Naegeli) Cleve-Euler. Besides that, 
in sites S0, Og and B1, oligotrophic and oligosaprobic 
taxa were identified such as Achnanthidium  jackii 
Rabenh and Surirella helvetica, with more eutrophic and 
beta-saprobic species like Amphora pediculus (Kütz.) 
Grunow (Van Dam, Mertens, and Sinkeldam 1994). At 
site In, meso-eutrophic and  oligo-saprobic taxa, 
Encyonopsis microcephala (Grunow) Krammer and more 
beta-saprobic and eutrophic taxa like Gomphonema 
minutum (Ag.) Agardh f. minutum (Van Dam, Mertens, 
and Sinkeldam 1994) were common. S1 and B2 sites 
had hypereutrophic and polysaprobic taxa such as 
Nitzschia palea, and Nitzschia capitellata Hust. (Van 
Dam, Mertens, and Sinkeldam 1994). Brackish and 
eutraphentic species were dominant in sites S3 and S6, 
as Craticula halophila (Grunow ex Van Heurck) Mann 
and Navicula erifuga Lange-Bert. (Lange-Bertalot 2001). 
Moreover, Nitzschia filiformis, known to be a species that 
can support brackish waters and tolerant of strongly 
polluted conditions, but Cyclotella meneghiniana Kütz. 
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Table-5: Indicator species analysis presenting only those species having results greater than 80% in a 
category. A. Relative abundance, which is the percentage of each species in each group defined by 
topographic position. B. Relative frequency, which is the percentage of samples in each group in which a 
species is present. C. Indicator values (% of perfect indication) which are a product of the relative 
abundance and constancy of a species in each group. The species marked in gray are found in A. B. and C. 


A. Group 1 Group 2 Group 3 
Achnanthes coarcta 100 0 0 
Fragilaria capucina 100 0 0 
Gomphonema clavatum 100 0 0 
Gomphonema angustatum 13 0 
Surirella biseriata 13 0 
Gomphonema olivaceoides 14 0 
Gomphonema gracile 0 20 
Nitzschia umbonata 88 12 
Hantzschia amphioxys 87 7 
Nitzschia scalpelliformis 0 89 
Navicula lanceolata 0 96 
Nitzschia recta 0 93 
Nirzschia sigma 0 95 
Navicula tripunctata 0 96 
Nitzschia obtusa 0 100 
Diploneis oblongella 0 100 

B. 
Gomphonema olivaceoides 100 50 0 
Navicula cryptotenella 100 50 0 
Gomphonema parvulum 100 100 0 
Amphora pediculus 100 100 67 
Hantzschia amphioxys 100 33 
Nitzschia fonticola 100 67 
Nitzschia dubia 100 67 
Craticula accomoda 100 67 
Cymbella affinis 100 0 100 
Mayamaea atomus 0 100 
N. recta 0 100 
N. sigma 0 100 
Navicula tripunctata 0 100 
Actinocyclus normanii 50 100 
Navicula veneta 50 100 
Cyclotella meneghiniana 50 100 
Diploneis oblongella 0 100 

C. 
G. olivaceoides 7 0 
Hantzschia amphioxys 87 2 
Nitzschia recta 0 93 
Nitzschia sigma 0 95 
Navicula tripunctata 0 96 
Diploneis oblongella 0 100 


and Navicula gregaria which support some pollution but 
not critical levels of pollution (Taylor, Harding, and 
Archibald 2007) were also collected. 


Conclusion 


The distribution of diatom species along 
environmental gradients has been the focus of numerous 
studies, which have been used to determine the structure 
of communities (Wellborn, Skelly, and Werner 1996; 
McPeek and Brown 2000). It has underlined that diatom 
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assemblages are driven by local environmental variables 
rather than regional processes (Soininen 2004). It may 
be useful, therefore, to describe such patterns and 
identify environmental drivers of species assemblages 
within different geographical settings, especially arid 
regions like North Africa. 

The present study attempts to identify the spatial 
distribution of benthic diatom assemblages and their 
responses to environmental gradients in a Moroccan 
hydro-system of the Sebou Basin. 
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The species distribution was mainly linked to human 
impacts and tidal influences in the estuary. The spring 
diatom communities of Sebou Basin are generally rich 
and diversified but numerically poor in distribution 
between taxonomic levels. Nevertheless, the distribution 
depends on the chemical quality of the water. Indeed, the 
stations downstream of anthropogenic activities have 
degraded water quality and low current flow aggravates 
the situation downstream from the city of Fez, which 
caused a decrease in taxonomic richness, diversity and 
stability of diatom community, dominated by pollution- 
resistant species, instead of pollution-sensitive species 
well represented in other unpolluted stations. It is hoped 
that the large treatment plant recently opened 
downstream of Fez will be beneficial to aquatic life of the 
Sebou Basin. 
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